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Abstract: In D = 4, A/" = 2 theories on M 3,1 , the index receives contributions not only from 
single-particle BPS states, counted by the BPS indices, but also from multi-particle states made 
of BPS constituents. In a recent work [1], a general formula expressing the index in terms of the 
BPS indices was proposed, which is smooth across walls of marginal stability and reproduces 
the expected single-particle contributions. In this note, I analyze the two-particle contributions 
predicted by this formula, and show agreement with the spectral asymmetry of the continuum 
of scattering states in the supersymmetric quantum mechanics of two non-relativistic, mutually 
non-local dyons. This provides a physical justification for the error function profile used in the 
mathematics literature on indefinite theta series, and in the physics literature on black hole 
partition functions. 
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1. Introduction 

The recent work [ 1 ] proposed a general formula for the index I(R, u , C) in four-dimensional 
field theories on M 3,1 with J\f = 2 supersymmetry. This index can be understood as the parti¬ 
tion function on M 3 times an Euclidean circle of radius R , with periodic boundary conditions 
for fermions, chemical potentials C conjugate to the electromagnetic charge 7 , and with an 
insertion of a suitable four-fermion vertex so as to saturate fermionic zero-modes. Equivalently, 
it can be defined as a trace 

T(R,u,C) = -|Tr w(u) (—1) 2 J 3 (2 J 3 ) 2 a, e~ 2 * RH - 2 ^’ c \ (1.1) 

over the full Hilbert space R[u) of the theory on M 3 (here J 3 is the angular momentum operator 
around a fixed axis, and cr 7 is a charge-dependent sign, satisfying the quadratic refinement 
property cr 7 oy = (—l)^ 7 ’ 7 , V 7+7 / ). 

Unlike the BPS indices 12 ( 7 ,u), which count single-particle BPS states and exhibit dis¬ 
continuities across walls of marginal stability, X(R,u, C ) is a smooth function of the Coulomb 
branch moduli u, away from the loci where additional states become massless. This is possi¬ 
ble because I(R, u, C) receives contributions not only from single-particle BPS states, but also 
from the continuum of multi-particle states. Indeed, while multi-particle states do not saturate 
the BPS bound M = |Z 7 |, the density of bosonic and fermionic states are not necessarily equal, 
as noted early on in [2, 3] (see [4, 5, 6 ] for recent discussions in the context of two-dimensional 
superconformal field theories). Still, only multi-particle states made of BPS constituents can 
contribute, so one expects that the index can be expressed in terms of the BPS indices f 2 ( 7 , u). 

In [1] , using insight from the study of the hyperkahler metric on the Coulomb branch in the 
theory compactified down to three dimensions [7], and by analogy with a similar construction 
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in the context of the hypermultiplet moduli space in string vacua [ 8 , 9], we proposed a general 
formula for the index 1 


R 


R 


d t 


I(R,u,C) = -—Z(X A F A ) + — Y, Sl (~t’ u ) I T (f-% - tZ,) log(1-*.,(()). (1.2) 


where £ 7 are the BPS rays {£ G C x : Z 7 /i G iM } , and T 7 are the solutions to the system of 
integral equations [7] 


X 1 = X^ exp 




47ri 


‘ d t't + t' 

YT 


7 log (1 - Xy(t')) 


(1.3) 


with X* { providing the ‘semi-flat’, large R approximation to X~ n 

_ a ^ e -7rm(t- 1 Z 7 -tZ 7 )-27ri(7,C) ^ ^ 

The fh 7 ’s are holomorphic functions on the twistor space Z of the Coulomb branch 
which provide canonical Darboux coordinates for the holomorphic symplectic structure on Z. 
They can also be understood as vevs of certain infrared line operators [11]. In the limit R —> oo, 
a formal solution to the system (1.3) is obtained by substituting X 1 —>- X* { on the r.h.s. and 
iterating. This leads an expansion of the form 


fh 7 = X^ exp 


Yl (“i.Q'i) n Cl(ai,u)g T 

T (ij)eTi ieTo 


(1.5) 


where T runs over trees decorated by charges cq such that 7 = 9t are certain iterated 

contour integrals [7, 12], and 0(7 ,u) = J0 d | 7 ^ 2 ^( 7 / d,u) are the ‘rational BPS indices’ [13, 14, 
15], which arise from expanding the log in (1.3). Substituting in (1.2), one obtains a formal 
expansion 

X = j(°) + X« + 1®' + • • • (1.6) 

7 7 > 7 ' 

where X (0 ^ stands for the first term in (1.2), while X 7 ").. )7ra , proportional to fl( 7 i, u)... 0 ( 7 n , u) 
is interpreted as the contribution of a multi-particle state of charge { 71 ,... 7 n } to the index. 
In particular, the one-particle contribution is obtained by replacing At, —y AX f in (1.2), leading 
to 

X« = ^a 7 n( 7 ,u) |Z 7 |iP 1 (2 7 ri ? |Z 7 |)e- 2 ^> (1.7) 

In [1], we matched this result with the index of a relativistic particle of charge 7 and mass 
|Z 7 |. To define the index, we regulated the infrared divergences by switching on a chemical 
potential 6 for the rotations J 3 in the xy plane, restricting the z direction to a finite interval 
of length L , and removing the regulators as follows, 



2 R lim d% 

0 -> 2 tT 
L—y 00 


"sin 2 ( 0 / 2 ) 
7tL 


Tr Ml( „, (cr, 


( 1 . 8 ) 


where R\(u) is the one-particle Hilbert space. The same regulator should then be used to define 
the full index (1.1). Our aim in this note is to perform a similar check for the two-particle 
contribution. _ 

1 Various refinements of this index have been introduced in [10], but lie beyond the scope of this note. 
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2. Two-particle contribution to the index 


According to the conjecture (1.2), the contribution of a two-particle state with charges { 7 , 7 '} 
to the index is obtained by inserting the one-particle approximation to (1.3) in (1.2), 


' 7 "( 2 ) _ 

X 7.7' ~ 


R 

647r 3 


(7,7')° ww) 


dt 

t 


7,7 


W 

7 


dt 't +1 ' 
t' t — t 


(t~% 


tz y ) <(«)*?(<') + (7 «■ 7') 

( 2 . 1 ) 


Defining -0 7 = argZ 7 , ip 77 i = ip 7 — t/y and changing the integration variables to t = ie 1 ^ 7 e x++x_ , 
t' = ie 1 ^ 7 'e x+_x_ , this can be rewritten as 


= ^(-D^ ( 7 , V) ^( 7 ) ^(Y) ^7+7' J dx + J dx - coth + ^77 

— OO — OO 

x [|Z 7 | cosh(x + + xS] + |Z y | cosh(i + - z_)] e - 2 »R|z, 7 o!ii(«++*-)- 2 *«|Zyl™h( I+ - I _)- 2 »i( 7+ y,C) i 

( 2 . 2 ) 


( 2 ) 

We shall be interested in the behavior of X) , in the vicinity of the wall of marginal stability 
where t/y 7 / —> 0, in the limit R —> 00 . We assume that 7 and 7 ' are primitive vectors generating 
the positive cone of BPS states whose central charges align at the wall, so that D( 7,77) and 
Cl(y,u) are constant across the wall, and equal to ^( 7 ,u) and ^( 7 ',«). Away from the wall, 
the integrals over x + and x_ are dominated by saddle points at x + = = 0 , producing 2 


' T ( 2 ) 

7 > 7 ' 


32vr 3 


- 1 ) (7,V> (7,7') ^7+7' 


IX 7/1 


l-Xy/1 


cot 



e -27rR(|Z 7 |+|Z 7 ,|)-27ri( 7 +7 , ,C> 


(2.3) 

However, the saddle point approximation breaks down near the wall where i\) 77 ' —> 0, as the 
saddle point collides with the pole at x _ = — 1 ' 0 77 '- To deal with this, we first perform the 
integral over x + , which is dominated by a saddle point at 


x + ~ 


| Zryl 

| - 

\Z 7 \ 

Z7I 

A | 

Z 7 ' | 


X- + 0 (x 2 _) . 


In the limit R —> 00 , X; 2 ,!, 


is well approximated by 


(2.4) 


x! 2) , 


(_1)(7,7'> (7,7') 0(7) 0( 7 ') cr 7+ y a/ i2(|Z 7 | + |Zy |) 


7,7 lg7r 3 


x / dx_ coth ( x_ + - ^ 77 / ) e 


4 irR\Z y \\Z , I (2 M 

-27ri?(|Z 7 | + |Z 7 ,|)- | Z7 ,;, 7, x 2 _-27ri(7+7',C) ^ ^ 


In the limit -0 77 / ->• 0, we can further approximate coth (a:) ~ 1/x and evaluate the integral 
using the formula [17, 4.18], valid for a and f3 real and non-zero, 


dz 


z — 1 a 


-e P 2 * 1 — i 7 rsgn (a) e a 2 /? 2 Erfc(| a/3 |) 


2 The analysis in this section bears some similarities with the one in [16]. 


( 2 . 6 ) 
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where Erfc is the complementary error function. Noting that 

1 


Zry | “ 1 “ jZ^'j r ^ J ^ Tfl^y ^ ^y , 

where m 7 y = is the reduced mass of the two-particle system, we find 

t( 2 ) ^ 


(2.7) 


~ I i- i (-l)<^'>n( 7 )f!(7')(7,7')<T 7 ^yfl(|Z 7 | + |Zy|) 
xsgnf^y) Erfc (\ip T /l 0rflm 77 .) e I- 2 ^i<y+y',C) 


( 2 . 8 ) 


The two-particle contribution is discontinuous across the wall: as 0 7 y goes from negative to 
positive, jumps by 

ai 7 v ~ T(-i)(^>(j( T )a ( y)( % y) ff7+l ,y fl (|z 7 | + |z 7 ,|)e- 2 " i ' z wi-2*‘<7+y',c) (2 9) 

On the other hand, the one-particle contribution Z] 1 ), is also discontinuous across the wall, 
due to the fact that the one-particle index 0(7 + 7') jumps [18] : 3 

0(7 + 7',w) = 0 + (7 + 7') - (-1) <7 ’ 7,> |(7,7') | 0(7)0(y) ©((7,7') ^77') • ( 2 - 10 ) 


The hrst term in (2.10) corresponds to the one-particle index on the side where (7,7') 0 7 y < 0 ; 
so that the two states of charge 7 and 7' cannot form a BPS bound state, while the second 
term is the contribution of the BPS bound state which exists on the side where (7,7') 0 7 y > 0. 
Inserting (2.10) in zj+y, and taking the limit R —* 00 , we find 


7"( l) 

“S+7' 


0 + (7 + 7 ') - (- 1 ) <7 ’ V> |( 7 ) l ') I 0(7) O(7 / )0((7,7') '077') 




\[W- 


87T 2 


2±Zl p -2 7 rfl|Z 7+7 ,|-2 7 ri(7+7 , ,C> 


( 2 , 11 ) 


whose jump exactly compensates (2.9). In fact, neglecting the difference between |Z 7+ y| and 
|Z 7 | + \Zy\ under the square root (which cannot be told apart in our approximation), the sum 
of (2.8) and (2.11) can be written as 


X (1 l , + Z (2) , 


|fl + (7 + 7') 

- ^(- 1 ) <7,7,> |(7,7')|0(7)n(y) ^1 + Erf (sgn ((7,7')) 0 7 y y/ 71R "©7')) } 


X (Tsy-^-ryf 


VW- 


7+7'I 


8tt 2 


e - 27 rR|Z 7 + 7 /|- 27 ri( 7 + 7 ',C> 


( 2 . 12 ) 


where we have used the identity Erf(x) = sgn (x) (1 — Erfc(|x|). In plain words, the addition 
of the two-particle contribution to Z^ , has converted the step function @(x) in (2.10) into 
the smooth function |[1 + Erf(x)]. This shows that the sum of the one and two-particle 
contributions is not only continuous, but also differentiable across the wall (see Figure 1 for 
illustration), which acquires a finite width of order 1/Rm lt y as a function of the relative 
phase ” 077 ' between the central charges Z 1 and Zy. It would be interesting to generalize this 
computation to the case of non-primitive wall-crossing, and to relax the non-relativistic limit 
R —»• oo. _ 

3 Here 0(x) denotes the Heaviside step function, equal to 1 when x > 0 and 0 otherwise. 
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Discretum Continuum Sum 



Figure 1: Behavior of the one-particle contribution (—20 (x)), two-particle contribution 
(sign(x) Erfc(|x|)) and their sum (—1 — Erf(x)) to the index with total charge 7 + 7 ' across a wall 
where the phases of Z 7 and Zy align (x -7 0). 


3. Supersymmetric electron-monopole quantum mechanics 


Our goal in the remainder of this note is to derive the two-particle contribution (2.8) from the 
supersymmetric quantum mechanics of a system of two non-relativistic particles with mutually 
non-local primitive charges 7,7'. After factoring out the center of mass degrees of freedom, 
which can be treated as in ( 1 . 8 ), and the internal degrees of freedom, counted by 
the system is described by Af = 4 quantum mechanics with Hamiltonian [19, 20 ] 4 

H = L (p ~ qA? ~ £ s ' d ' ® (l2 " aa) + L (*";)' ■ (31) 

where q = 1 ( 7 , 7 ') is half the Dirac-Schwinger-Zwanziger product of the electromagnetic 
charges, B = is the magnetic field of a unit charge magnetic monopole sitting at the 
origin, A is the corresponding gauge potential, a are the Pauli matrices, and m = rri~ n < is the 
reduced mass of the two-particle system. Classically, the system has bound states for qd > 0, 
no bound states for qi3 < 0, and a continuum of scattering states with energy E > E c = y± 
The parameter is fixed by equating E c with the binding energy, 

— = |Z 7 | + |Zy| - |Z 7+ y|, ( 3 . 2 ) 


so $ ~ m , 077 ' near the wall, cf. (2.7). Quantum mechanically, El describes two bosonic 
degrees of freedom with helicity h = 0, and one fermionic doublet with helicity h = ±1/2 and 
gyromagnetic ratio g = 4. This unusual value is fixed by the requirement of supersymmetry, 
and can be understood as the combined effect of electromagnetic and scalar interactions [22], 
Indeed, the Hamiltonian (3.1) commutes with the four supercharges (here II = p — qA ) [19, 
20 , 21] 


Q a 


Qa 


1 ( 0 -i - 2 ) + . n 

V2m \i (tf - l) ± cr • fl 0 


1 / 0 - (# - l) a a - in a ± e abc Il b a c 

a / 2 m ($ - r) a “ + in « + e abcA b a c 0 


(3.3) 

(3.4) 


4 A f = 4 supersymmetry allows a position-dependent rescaling of the kinetic term [21], but the spectral 
asymmetry is independent of this deformation, as long as it goes to one at spatial infinity. 
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which satisfy the algebra (here m — 1,2, 3,4) 

{Q m ,Qn} = ‘2H8 mn . (3.5) 

The complete spectrum of this Hamiltonian was analyzed in [23], but unfortunately these 
authors stopped short of computing the density of states in the continuum. We shall revisit 
this computation, adapting the classic treatment of the electron-monopole system without 
potential in [24]. 

The Hamiltonian (3.1) commutes with the total angular momentum operator [24] 

_ -> T 1 

J = f A (p- qA) - q- + -<J<g> (1 2 - <r 3 ), [J a , J h ] = kabcJc- (3.6) 

r 4 

The Schrodinger equation 11 + = E+ can be solved by separating the angular and radial 
dependence. For this we diagonalize J 3 and J 2 denoting by m and j (j + 1) their eigenvalues. 
For the spin 0 part (corresponding to the first two entries of the eigenvector T), we write 


* = /M y,j , m , j> |?|, j-q ez, 


(3.7) 


where Y q ^ m are the monopole harmonics (also known as spin-weighted spherical harmonics), 
given in the patch around 6 = 0 by [25] 


y _ 9 m 

1 q,l,m ^ 


(21 + 1) {l m)!(/ + m)! (i _ cos + cos ^ P ~™~ m ( cos 6) e< m+ ^ . 


4vr (l - g)! (I + q)\ 

Here P^ ,/3 (x) are the Legendre polynomials 


Using 


i— 1 i Tl d” 

P+(x) = ^-(1 - x)-“(l + x+— [(1 - *)-(! + xf**} 


_ 1 1 / ^ \ 2 
(p — qA) 2 = — d 2 r + — [r A (p - qA) 

ry* ryZ y 

\ 2 


r A (p — qA)j = L~ — q~ = 


-do sin Ode 


sin 6 u sin 2 6 

one can show that the radial wave function satisfies [25] 


(3.8) 

(3.9) 

(3.10) 

(9f + \q(cos6 - l)) 2 , (3.11) 


19 ,v‘-q 

- d r H- 

2m r r 2mr 2 


2 1 




2 m 


rJ 


f(r) =0, 


(3.12) 


where v — j + 

For the spin 1/2 part (corresponding to the last two entries of the eigenvector T), the 
angular dependence is a linear combination of modes with orbital momentum j — \ and j + 4 

[24], 


6 (1) 


j+m y i \ 

94 1 4 


2 i 


1 


q,j- 2 ’ m ~ 2 
J —Y 1 1 

2 1 q,j~ 2> m +2/ 


b® 

J,m 


(_ j j-m+l y , , \ 


2j+2 * q,j+\,m-\ 


j+m+l y 

\ V 2 t+ 2 q,j+ 2- m +2 J 


, j — q G Z + - (3.13) 
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The first set of modes occurs for j > \q\ + | while the second occurs for j > \q\ — In order 
to diagonalize the action of a ■ ■r'and a ■ (p — qA), which commute with J, it is convenient to 
introduce the linear combinations (for j > |g| + |) [24] 


t (+) _ r A 1 ) 

a> j,m t 7, 


c d> {2) 

j,m ' J,m ’ 


= C - <t> { Sn + C + 


3,™ 


3,m 


(3.14) 


where the coefficients 


c± = 


q (a /2 j + 1 + 2 q ± a /2 j + 1 — 2q) 


satisfy ci + c 2 _ = 1. Using r ■ a = 


191^2(4^ + 2) 

do,i,o \/2do,i,-i 

^ V-V^do,!,! — do,i,o 


(3.15) 


, and the multiplication rule 


ii+i2 


y . y . 

J <n 1 92 j2,m2 


E 

33 =max( I ji - J 2 1, | m -1 +"> 2 1) 


4vt(2j 3 + 1) 
x {ji,m 1 ,j 2 ,m 2 \j 3 ,m 1 +m 2 )Y qi+q2 ,33,mi+m2 


(2ji + 1)(2 j 2 + 1) , . . I . v 

{ji,-qi,j2,-q2\j3,-qi - 92 ) 


(3.16) 


for monopole harmonics, one can show that these combinations satisfy, for any /(r) 


a 


(a ■ r) /(r)4t! = ~r fir )^J, 

• (p~ qA) /(r)4t! =i (<^ + r_1 (l T /i)) /(r) 


where we defined 


p=\lij + \) 2 -q 2 - 


Using the fact that the Hamiltonian in the spin 1/2 sector can be written as 


H v* = i - Y)) 2 -£S.g + ±(t-i)‘ 


and the identity 


d r + 


1 ± /i 


r 

4 ±) : 


d r + 


1^) = -air 

r / r 


/i(/i y 1) 


we find that its action on /(r)y m is given by 
#1/2 • /(r) 


’£ (+) \ 

I _ 


e<-> i 

V s .?,"V 


—-<9 2 r + + — (■& - -Y 

2m r r 2 mr 2 2m \ r) 


f( r ) 


X 


d+) N 

sum 

(-) 


-/i q 


f( r ) 


x i+r 

A-) 

^>7,1 




(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


2 mr 2 y q /y 

The 2x2 matrix on the second line has eigenvalues ±y/u 2 + q 2 — ±(j T |), and eigenvectors 

f$ = (M =F #+?) 44 - 9 44 U 22 ) 
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Noting that the coefficient of the centrifugal 1/r 2 term in the potential is proportional to 


h- 2 ± vV + g 2 = (i + \ ± ^ - <? 2 ~ \ . (3-23) 

we find that the radial wavefunctions f± ( r ) associated to the eigenmodes satisfy the same 
equation as (3.12) with v — j +1 (for the + sign, which we refer to as the helicity h = | mode) 
or v = j (for the — sign, which we refer to as the helicity — | mode). 

Finally, for j = \q\ — 1/2, the space of eigenmodes of J 2 , J 3 is one-dimensional, spanned by 


j+m-l ( s i n # \ 

rim = 4 m (i + cos 2 i 1 - cos #) 2 e l(m+j) ^ _ cog ^ J • ( 3 . 24 ) 

One has, in place of (3.17), 

Q. Q. 

(a ■ r)rim = r—rim, a ■ (p - qA) f{r) rj m = -i—r (d r + r -1 ) f(r) q m , (3.25) 

lei lei 

leading to the same equation (3.12) with u = j. 

In summary, the radial equation is given by (3.12) with 

v = j + h-\-\, j = \q\+h + £ (3.26) 

with h — 0 for the two bosonic degrees of freedom and h = for the spin 1/2 degree of 
freedom, and £ e N in all cases. Solutions to (3.12) with energy E = ^ ^ are linear 

combinations 5 


rf(r) = (3 M m v (2irVk 2 — d'A + y W i q o (2i ry/W--£p) , (3.27) 

yj i?2 ’ V / fc2_$2 ’ V / 


where Mx, u (z) and Wy^z) are Whittaker functions, which are solutions of the second order 
differential equation 






u/y) = 0 


(3.28) 


satisfying 

M XtV (z) ~ ^ , Wy„(z) ~ yV z/2 . (3.29) 

-2—>-0 |,z|—^OO 

The solution proportional to M is regular at r = 0, while the solution proportional to W 
describes an outgoing spherical wave. Since 


MxM = W-Uj'z) + W ^). (3-30) 

r ( u - A + |) i(u + A+|) 


5 It helps to note that (3.12) is isomorphic to the Schrodinger equation of the hydrogen atom, whose radial 
wave-functions are linear combinations of M. i (2ifcr) and W. i (2ifcr) where q = qxqivnlk , E = k 2 /2m. 

i 7 j 2 1??,2 

(see e.g. [26, Chap. 14.6]). 
















we find that the S-matrix in angular momentum channel i and hclicity h is 


s h M = 


V(v + A + \) _ r (M + t + 2h + 1 - iyj=f) 


r(^-A+ 2 ) Y (\q\+£ + 2h + l + i- 


q-d 


(3.31) 


Vfc 2 -!? 2 , 

In particular, bound states correspond to poles of the S-matrix, and occur only when q& > 0, 
with energy 


■d 2 

Eh,i,n — - I 1 — 


q- 


2m \~ (\q\ + i + 2h + n + l) 2 ^ 

with n > 0. The energy depends only on the sum N = £ + n + 2 h, so the spectrum has 
additional degeneracies beyond those predicted by rotational symmetry and supersymmetry 
[19, 23]. The supersymmetric ground state occurs in the h — — | sector with n = l = 0 and 


has degeneracy 2\q\ = 1 ( 71 , 72 )!. Its wave function is 


-q—1 p — fir. 


V m 


, in agreement with [20, 4.16]. 


The density of states (minus the density of states for a free particle in M 3 ) is the derivative 
of the scattering phase, p[k) d k = ^d log S(k). The canonical partition function for states of 
helicity h is therefore 


Tr h e 


-27 tRH 


e M) + %| + 2 h + 1) e~ 2irREhA 

0 d kd k 


i=0 n =0 


^ ((21 + 2\q\ + 2 h + 1) 


e=o 


'H 


27ri 


log 


T (jg| + £ + 2h + 1 -i^f=) 


T (|q , | + £ + 2 h + 1 + i 




y/k 2 — 19 2 


(3.33) 


where the first term, corresponding to discrete bound states, contributes only when gd > 0. 
Summing over all types h weighted by fermionic parity (—l) 2/l , only the BPS bound state with 
E — 0 contributes from the first term, while the contribution of the continuum of scattering 
states simplifies to 



(2t + 2|g|) log — 

z t 


(2£ + 2|g| + 2)log^±i 

z i +1 


e 


nRk 2 

m 


(3.34) 


where 


z e = \q\ + £ - i- 


q£} 


Vk 2 - E 2 ‘ 

Cancelling the terms in the sum, only the contribution l = 0 remains, leading to 

I 91 h/t.2_,<)2 J nRk 2 2 qi£ 


m 


dkd k 

~2^r log 


k=[ 


\Q\ +i 


\/k 2 — 1? 2 

<g 

Vk 2 -tf 2 


7T 


fc=h 


dfc 

k^/k 2 -® 2 


(3.35) 


(3.36) 


This is in fact the standard spectral asymmetry predicted by Callias’ theorem [27] 6 . In¬ 
deed, the result (3.36) does not depend on the details of the S-matrix, but only on the ra- 

6 While the Callias theorem is valid a smooth monopole background, the extension to singular monopoles 
was worked out in [28], and leads to the same spectral asymmetry. I thank A. Royston and D. van den Bleeken 
for discussions on this matter. 
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tio Sh/+i(k)/Sh/{k), which as we discuss in an Appendix, is fixed by supersymmetry in the 
asymptotic region. 

Using 


OO 



fc=|#| 


d k _ ■rrRk ^ 

-(3 m 

kVT? - tf 2 


2|i?| 


Erfc 




(3.37) 


and adding in the bound state contribution, we get, for arbitrary signs of q and t), 


Tr(—l) F e wtH = — \2q\ Q(q'd) + sign(g'd) |g| Erfc 


= -q 


sign(g) + Erf 




(3.38) 


This is indeed a smooth function of d, which interpolates from 0 at'd = — oo to —2 q at i9 = +oo 
when q > 0, or from —2\q\ at = — oo to 0 at t? = +oo when q < 0 (see Figure 1, which displays 
the case q = 1). Including the degeneracy 0 ( 7 ) 12 ( 7 ') of the internal degrees of freedom, and 
the contribution of the center of motion degrees of freedom, given in the last line of (2.12), we 
find perfect agreement with the two-particle contribution to the index predicted by the formula 
( 1 . 2 ). 


4. Discussion 

In this note, I have shown that the general formula for the index (1.2) in J\f = 2, D = 4 
gauge theories correctly reproduces the contribution of the continuum of two-particle states, 
in the vicinity of a wall of marginal stability where the constituents can be treated as non- 
relativistic BPS particles. In particular, I demonstrated that the contributions of the BPS 
bound states and of the two-particle continuum add up to a smooth function of the moduli, 
even though each of them is separately discontinuous across the wall. This analysis provides 
a physical justification for the replacement sgn (x) —>■ Erf(x), which has been postulated in 
studies of black hole partition functions in order to enforce S-duality or modular invariance 
[29, 30, 31, 32], a trick borrowed from the mathematics literature on indefinite theta series 
[33]. It would be very interesting to calculate the contribution of the continuum of multi¬ 
particle states away from the wall, a challenge that will require to understand the dynamics of 
a collection of relativistic mutually non-local particles beyond the BPS regime. 

It is worthwhile noting that similar smooth interpolations across walls of marginal stability 
have been encountered recently in localization computations of the index in gauged supersym¬ 
metric quantum mechanics in certain scaling limits [34, 35]. More generally, error function 
profiles are ubiquitous in the context of Stokes phenomenon [36], which is formally similar 
with wall-crossing [37]. It would be interesting to explore these connections. 

Acknowledgments: It is a pleasure to thank A. Neitzke, J. Manschot, G. Moore, A. Royston, 
A. Sen, J. Troost, D. van den Bleeken and P. Yi for useful discussions. Special thanks are due 
to S. Alexandrov for collaboration at an initial stage of this work and continued collaboration 
on closely related topics. 
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A. Robustness of the spectral asymmetry 


In order to elucidate the origin of the cancellations in (3.34), we need to understand how 
supersymmetry relates the density of states in the bosonic and fermionic sectors. For this 
purpose, notice that the operators 


Qr =d z - 

Q'r = 9 z + 


A 

2/7 + 1 

A 

2/7- 1 


1 

2 i Qk 2 - d 2 
1 

2iQk 2 - d 2 


d r — 
d r + 


v 


v ~\ 


+ 


q& 

v +\. 

qd 


(A.l) 


maps solutions of the Whittaker equation (3.28) with parameters (A, v) to solutions of the 
same equation with parameters (A, v + 1) and (A, v — 1), 


Qr ■ T^\,v ~ T>X,v+l ■ Qr , Q' r ' B>x,v ~ T^X.u-l ' Q' r 


(A.2) 


In fact, Q r and Q' r can be interpreted as the supercharge for the radial problem. To see this, 
consider acting with y/2m Q1/2 = <7 • fl + i(® — 1 ?) on fermionic eigenfunctions f± '■ this 
should produce linear combinations of bosonic eigenfunctions with the same energy and spin 
j ± |, namely f\ cp'^ m and 02 dfln- H er + /+, /-, fi, $2 are solutions of the radial equation (3.12) 
with v — j + 1, j, j, j + 1, respectively. Indeed, we find 


\\Phn Q a • /+ = -( 2 j + 1) c- (d r + t) /+ 4$L + 2 ^ c - f+ 

iV2m Q 4 • /_ = -(2 j + 1) c+ (d r - J —^ r 1 + _ 2 /^ c + /_ 0^ 


(A.3) 


The differential operator in brackets coincides with r" 1 ■ Q' r ■ r and r' 1 ■ Q r ■ r, up to overall 
normalization. Acting on the Whittaker wave-functions, we have, as a consequence of (A.2) 
and (3.29) 7 , 


Qr ' W X}U (Z) 
Qr • W.^(-Z) 

Q r ■ M x , v {z) 


j 

2/7+ 1 
A + Z 7 + A 

2 // + 1 


WW*) 

W_ A> , +1 (-+) 


(^ + |) 2 -A 2 

2(2/7 + i)2(/7 + i) 




(A.4) 


It follows from these relations that the reflection coefficients A(A, /z), B( A, //) defined by 


M\ iV {z) — A(A, v) W- X , v 


(-z) + B(\,v)j<- v -&W x , v (z) 


(A.5) 


satisfy 

A(A, v + 1) 2(/7 + 1)(2z7 + 1) B( A, /7 + 1) 2(/7 + 1)(2z7 + 1) 

A(A,/z) /7 - A + \ ’ B(\,v) z7 + A + ^ 


7 The last equation in (A.4) was noted in [23, VI.24], 


































Denoting the reflection coefficient S„( A) = A(X,u)/B(X,u), one has 


^+i(^) v + X + \ 

57(A) v-X+\ 1 ‘ j 

This relation only depends on the first two equations in (A.4), which in turn follow directly 
from the action of Q r ~ d r + on the leading asymptotic behavior W\ iV ~ z x e ~ z A of 
the incoming/outgoing plane waves. More generally, the ratio (A.7) depends only on the 
supercharge at radial infinity and should be unaffected by short-distance corrections to the 
potential or to the conformal factor in the metric. 
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